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Abstract 

Charge asymmetry in electron (positron) scattering arises from the interference of the Born 
amplitude and the box-type amplitude corresponding to two virtual photons exchange. It can be 
extracted from electron proton and positron proton scattering experiments, in the same kinematical 
conditions. Considering the virtual photon Compton scattering tensor, which contributes to the 
box-type amplitude, we separate proton and inelastic contributions in the intermediate state and 
parametrize the proton form-factors as the sum of a pure QED term and a strong interaction term. 
Arguments, based on analyticity, are given in favor of cancellation of contributions from proton 
strong interaction form factors and of inelastic intermediate states in the box type amplitudes. In 
frame of this model, and assuming a dipole character of form-factors, numerical estimations are 
given for moderately high energies. 

PACS numbers: 
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I. INTRODUCTION 



annihilation channels 
proton form factors (FFs) 



Recently, a lot of attention was devoted to 27 exchange amplitude both in scattering and 

2 1, in connection with the experimental data on electromagnetic 

Extraction of box type (two photon exchange amplitude (TPE)) contribution to elas- 
tic electron-proton scattering amplitude is one of long standing problems of experimental 
physics. It can be obtained from electron proton and positron proton scattering at the 
same kinematical conditions. A similar information about TPE amplitude in the annihila- 
tion channel can be obtained from the measurement of the forward-backward asymmetry in 
proton-antiproton production in electron-positron annihilation (and the reversal process). 

The theoretical description of TPE amplitude is strongly model dependent. Two reasons 
should be mentioned: the experimental knowledge of nucleon FFs is restricted in a small 
kinematical region and the precision of the data is often insufficient, and the contribution 
of the intermediate hadronic states can be only calculated with large uncertainty. 

A general approximation for proton electromagnetic form-factors follows the dipole ap- 
proximation: 

G E (q 2 ) = = G D (Q 2 ) = (1 + Q 2 /0.71 GeV 2 )- 2 , Q 2 = -q 2 = -t, (1) 

where /1 is the anomalous magnetic moment of proton. However, recent experiments [4] 
showed a deviation of the proton electric FF from this prescription, when measured follow- 
ing the recoil polarization method, which is more precise than the traditional Rosenbluth 
separation. Such deviation was tentatively explained, advocating the presence of a two 
photon contribution. 

The motivation of this paper is to perform the calculation of charge odd correlation 

dd _ da^ - da^ 
A ~ 2dbf ' {) 

in the process of electron-proton scattering in frame of an analytical model (AM), free from 
uncertainties connected with inelastic hadronic state in intermediate state of the TPE am- 
plitude. In frame of this model it is possible to show that the effects due to strong interaction 
FFs and those due to the inelastic intermediate states almost completely compensate each 
other, within an accuracy discussed below. 
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Our paper is organized as follows. In part II, III we consider the charge-odd contribution 
of triangle and box-type diagram. In part IV we describe the procedure of numerical inte- 
gration. In In part V we present the results of numerical integration for asymmetries and in 
the Conclusions we estimate the accuracy of the obtained results. The appendices contain 
the tables of four-fold integrals and some details of calculations. 

II. ANALYTICAL MODEL FORMULATION 

In the analysis of the TPE amplitude we consider the electromagnetic interactions in 
the lowest order of perturbation theory. Hadron electromagnetic FFs are functions of one 
kinematical variable, Q 2 and the static value of the Dirac FF of the proton (for Q 2 =0) is 
unity due to QED origin. Therefore we parametrize the proton FFs in the form 

F^q 2 ) = 1 + F ls (q 2 ), F 2 (q 2 ) = F 2s (q 2 ), F ls (0) = 0, F 2s (0) = p, (3) 

Let us discuss now the arguments in favor of a cancellation of the terms of order of F 2 with 
the contribution of the inelastic hadronic intermediate states, in TPE amplitude. 

The TPE amplitude contains the virtual photon Compton scattering tensor. It can be 
splitted in two terms, when only strong interaction contributions to Compton amplitude 
are taken into account. One term (the elastic term) is the generalization of the Born term 
with the strong-interaction FFs at the vertexes of the interaction of the virtual photons with 
the hadron. We suppose that the hadron before and after the interaction with the photons 
remains unchanged. The second term (inelastic) corresponds to inelastic channels formed 
by pions and nucleons or similar hadronic states which can be excited in the intermediate 
state, between the vertexes of the virtual photon interaction with the hadron. 

For this aim let us present the loop momentum integration element in the form 

d 4 A; = ^-d 2 k ± d Sl ds 2 , (4) 

where s = 1pp\ is the total energy, s± — (k — pi) 2 , and s 2 = {k + p) 2 are the invariant mass 
squares of the upper (electronic) part of TPE Feynman diagram and the lower (hadronic) 
ones, d 2 A;j_ represents the integration on the components of the loop momentum k±pi = 
k±p = which are transversal to the initial electron pi and proton p momenta. 

We can consider the upper and the lower block tensors to be both gauge invariant (the 
factor 1/2 is introduced to avoid double counting and two Feynman diagrams are included in 
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(N) (N+7& (N+27C) 

(a) 



(2NN) 




each block). Therefore, the TPE amplitude can be written in the form (we omit the factors 
corresponding to fermion spinors): 

1 (47m) 2 (27n) 2 f d 2 k ± ds 2 

where L^ v = and is the Compton tensor of proton. The integration contour 

is drawn in Figfjji: it starts from — oo — iO and follows to +00 + iO, so that it belongs to the 



physical sheets of s and u channels [5]. 

In the physical sheet, the Compton amplitude has a pole, corresponding to a single proton 
state in the intermediate state and two cuts: the right one, corresponding to the inelastic 
states in S2-channel and the left one, starting at s 2 < — 9M 2 (M is the proton mass). Closing 
the integration contour to the left and to the right side (see Fig. [Tfci) it was shown (see Q]) 
that the following relation holds 

A-left -^-elastic -^inelastic' (6) 



Comparing with the sum rule obtained in [5], here the Born amplitude contribution is 
omitted, as well as only strong interaction effects are considered here. Omitting the left 
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(c) 

FIG. 2: Compton scattering Feynman diagram illustrating the left cut (a) which is equivalent to 
the u-channel discontinuity of the Feynman amplitude (b); The gauge invariant set of Compton 
subdiagrams is shown in (c). 

cut contribution Ai e f t , (our estimate shows that it can be included in 10% error bar j^]), 
the effects of strong interaction contributions to the hadron FFs compensate the strong 
interactions contributions arising from the inelastic channels. 

As an example, the Feynman amplitude at the origin of the left cut in the S2 plane of 
virtual Compton scattering which contributes to the TPE blocks, is drawn on Fig. [2^, un- 
derlying the proton propagators which corresponds to real 3-proton u-channel intermediate 
state. 

Its physical meaning is the interference of amplitudes of proton- antiproton pair pro- 
duction in virtual photon-proton collisions due to Fermi statistics (see Fig. [2Jd). A rough 
estimate of the ratio of contributions of typical right hand cut R T % g ht to the left hand cut Ruft 
is the ratio of cross sections of pion photoproduction to nucleon-antinucleon photoproduction 
cross section on proton: Rieft/ Rright ~ (2Mm 7r )/(10M 2 ) < 10%. 

This argument was shown to be exact in frame of QED |||, where a rather specific 
kinematics was considered: forward scattering amplitude at high energies. The application 
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to the case of nonforward TPE amplitude requires a more rigorous proof, which is outside 
the purpose of this paper. Here we suggest to consider our approach as a model, the validity 
of which should be experimentally verified. 

For example, experiments measuring charge-odd observables in ep scattering will be crit- 
ical for verification of the validity of our model. 

Proton FFs enter in the box amplitude in a form which can be schematically written as 

f d*k 1 + F s (k 2 ) 1 + F s (k 2 ) 

J Z7r 2 (e)(p) (A:) (*) ' U 

(A;) = k 2 - A 2 , (k) = (k - qf - A 2 , (e) = (k - Pl ) 2 - m 2 , (p) = (k + p) 2 - M 2 , 

where we extract the QED part and do not distinguish Dirac and Pauli form factors. This 
expression can be rearranged as 

1 + F s (k 2 ) 1 + F s (k 2 ) F s (k 2 )F s (k 2 ) 1 \F s {k 2 ) F s (q 2 



(k) 



(k) 



+ 



F s (k 2 )F s (k 2 ) J_ 
(*)(*) + (k) 
F s {k 2 ) F s {q 2 



(k) 



Q 



(k) 



(k) 



+ 



F s (q 2 



1 1 

W) + W) 



(k)(k) 



A Int + A FR + A TR + A 



l-Box- 



(8) 
(9) 



According to our model assumption the first term Ai nt in the right hand side (r.h.s.) of 
(Q is compensated by the inelastic intermediate hadron state contribution and so it will 
be omitted below. Next two terms in (jSj) do not contain infrared (IR) singularities but 
contain the ultraviolet (UV) ones. The fourth term in (jHJ) suffer from both IR and UV 
divergences, the last one suffers only from the IR divergences. Contributions of the last two 
terms can be calculated analytically as well as it does not contain FFs uncertainties at the 
loop momentum integration. The explicit results for them are given below. 

Keeping in mind the UV convergence of the initial amplitude, we extract the UV cut-off 
A depending contribution containing In from the fourth term in (jSJ) and add it to the 
contribution of the 2 and 3 terms providing their UV convergence. This procedure denoted 
in equation (J9j). Part of calculations concerning pure QED contribution (the last term in 
this equation) was performed in our paper [2|. 

The cross-section of elastic ep scattering 



e(pi) +p(p) e(pi) +p(p) 



(10) 
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in Born approximation in laboratory frame (p = (M, 0, 0, 0)) has a form: 



with 



a 2 cos 2 \ 1 



dtt e(l + r) 



2E . 2 e 



AE 2 sin 4 § p 



~j P — 1 + TV sin r 



£il — ^1, Q 2 = - g 2 = -t = 2pip j , s = 2ME, u 



\f "'2' AM 2 ' 
s 
P 



-2pp[, 



n 

s + t + u = 0, e' 1 = 1 + 2(1 + r)tan 2 -, 



11^ 



&red — t~G 2 m + eG 2 E , Gm — Fi + F 2 , Ge — Fx — tF 2 . 



(12) 



Here 6- electron scattering angle and F\ = 1 + -Pis, -F2 — -Pis - Dirac and Pauli proton's FFs. 

Keeping in mind the representation in the form of Eq. it will be convenient to use 
another (equivalent) form of Born cross section: 



da B da Bt da Bbox 



a 



dn dn 



dn 



M 2 pH- 



:{B t + B, 



box I 



B t = i(F 2 - F x )(2tM 2 + s 2 + u 2 ) + t 2 F 1 F 2 - rF 2 (tM 2 + su) - ^F 2 t 2 , 
B box = ]^F l {2tM 2 + s 2 + it 2 ) + h 2 F 2 . 



(13) 



(14) 



The IR divergence from virtual photon emission contribution is, as usually, canceled when 
summing with contribution from emission of soft real photons: 



da soft 
dQ 



da Bt da Bbox ] M da 

-dn + —~ l<) 



soft 



soft 



J s 

Bt , aa Bbox 



The quantity S°^ t was considered in 



dn \ soft 



dn 



dn 



(15) 



xodd 
°soft 



167T 3 J uj p[k pik p'k pk 



S ,u<AE 



2a 

7Y 



' 1 , 2pAE 
2 In - In -4 — 

P A 



In x In p + Li 2 ^1 j — Li 2 ^1 



x 



with A, AE- soft photon mass and it's maximal energy in laboratory frame. 

Note that d°^j t does not have a definite symmetry under substitution p <-> — p'\ s 
due to the specific definition of soft photon in the laboratory frame. 



(16) 



u 
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The virtual contribution to the cross section can be splitted in three terms (according 
Eq. ©): 



da v da vt da vb da F 



a" 



: (a t + a b + a f ). 



(17) 



dtt dtt dtt dtt 2nt 2 M 2 p 2 
The first term appears from the contribution of triangle-type diagram and can be put in the 
form: 



a t 



[l-P(s^u)) 



d A k 1 

in 2 (k) 



S e Sj> S S S T 



(e)(p) (e)(p) 



with 



(e) = (k + p\) 2 - ml (p) = (p' - k) - M 



(18) 



(19) 



where P(s <-> u)f(s,u) = f(u,s) is the substitution operation and 
1. 



S e 



Trp'^ip! - fc)7^i7r?, 



S T = -TrR{(T^q)- ltl )(p + k + M) lv 



Sn 



-TrRjJp' -k + M){Y v {q)- lv ) 



?M) = l^q 2 ) + m F 2 (q 2 )(q ltl - ^q); R = (p + M)T v (-q){p' + M). (20) 



The box type contribution is parameterized as: 



a b 



[1-P{8 
1 



«)) 



d 4 A; 



in 2 (jfc)(jfe)(p)( e ) 



(21) 



„ - -TrR 1 ^(p + k + M) lp 



Some of the necessary integrals have been previously calculated in [6]. However, for com- 
pleteness, they are all given in Appendix A. 

The finite part contribution (2nd and 3rd terms in rhs of Eq. (Q) are parametrized as: 



a f 



:i-p(s 



«)) 



d A k t 



m 2 (k) L(e)(p) (e)(p)J 



- A f (u,s), 



(22) 



and 



5 F 



TrR^(p + k + M)<f> v ; 
TrRQJp 1 -k+ M) lv , 



7m 



4M 



F 2 (fc 2 



.(23) 



In the integration over the four- vector k in a/ and at enters the UV cut-off parameter 
\k\ 2 < A 2 . To obtain Apr and Atr in Eq. (Q the terms containing ln(A 2 /M 2 ) are singled 
out from the contribution of the fourth term of rhs (jHD and add it to the contribution of 2nd 
and 3rd term in (jSJ). So we can put in (17) a t + a^ + af = atr + ctb + dfr, explicitly eliminating 
the cut-off parameter dependence. 

III. VIRTUAL AND SOFT PHOTON EMISSION CONTRIBUTIONS OF TRIAN- 
GLE TYPE DIAGRAM 

The interference of Born amplitude with the part of TPE arising from the fourth term 
from the r.h.s. of Eq. (jUJ) with the corresponding part of soft photon emission leads to 

da T da vt da S g t ft 



dn dn 




— In 2 p H — hip In x H — Li 2 I 1 



px 




(F 2 + F 1 )(l-F 1 )ln^, 



(24) 



with 



D 



tsv 



[, s o — u ( M 2 \ ( M 2 \ „ 




(25) 



and 



A(s, t) = aiFi(l - Fi) + a 2 F 2 + a 3 F| + a 4 FiF 2 , 



B(s, t) = 6iFi(l - F x ) + 6 2 F 2 + 6 3 F 2 2 + 6 4 FiF 2 , 
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with 



a i 

a 2 



+ ti ( 3M2 + 4s ) + 7h 6M4 + 4sM ' + 2tM ' + 6u "] 
2[M Z + s) z 4 



tM 4 (7M 2 + 9s) t 
4(M 2 + s) 2 4 
M 2 i 



(5w + 7M 2 



(2M 4 + 5st + 3sM 2 + 5tM 2 ) + — -(2M 4 + M 2 (w - 2t) - 8su 

16M Z 



16(M 2 + s) 2 

M 2 t t 

(8M 4 + 2st + 23sM 2 + 2tM 2 ) - -(-18M 2 + 18t + 13s), 



8(M 2 + s) 



»,2n M\M 2 +t) , 7 tM 4 



si 



+ 



tM 4 



15 5tM 4 
o 4 = — st + 



(26) 



16 16(M 2 + s)' "* 8"" ' 8(M 2 + s) 

Note that function D ts b in r.h.s. of Eq. (!24|) changes the sign at substitution s <-> w, in 
particular 

2 — S — zO 



S 1 D 

p = > -, He 

—u p 



M 2 



—u 



—u 



In 2 — -tt 2 ^ln 2 — ; In— ^ln— . (27) 
M 2 M 2 ' M 2 M 2 y ! 

The last term in (24), which contains the cut-off parameter is necessary in order to remove 
the A-dependence of a/. 



IV. VIRTUAL AND SOFT PHOTON EMISSION CONTRIBUTIONS OF QED 
BOX TYPE DIAGRAM 



The box-type contribution (last term in (jHJ)) with corresponding part of soft photon emis- 
sion is given by (the list of necessary loop momentum integrals and details of the calculation 
are given in Appendix A): 



da B 

an 



da vb + da s °£ x 



dQ 

2a 3 
TiM 2 p 2 t 2 

+ Li 2 | I 



dQ e 

B b ox 
1 

px 



p 4tM 2 



Li 2 1 



2 In p + In p In x 

- t \-P(s^u)](d 1 F 1 -d 2 F 2 ) 



(28) 
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with 



\ ln 2 (4r) - JL- ln(4r) + M 2 Fq Ut + 2 - J5£-) 



, 2 s 9 ( M 2 

In 2 _ + vr 2 + 2Li 2 1 + 

—t \ s 



l-2r) 



do 



4r 

s 



2 ^ 2 — - I - M 4 



21npln(4r) - In 2 — + 7r 2 + 2Li 2 1 + 

\ S 



2M 2 



su\ ln^fa 



2(1 + r) 
-Li a (1 + 



2M 2 



ln(4r) 



M 2 



M 2 + s 



M 2 



t J s + M 2 ' 
,„^-,„„„(4 T) + i(.n^-^) 

(29) 



and Fq given in (14T|) . 



V. INFRARED SINGULARITIES FREE CONTRIBUTIONS OF BOX AMPLI- 
TUDES 



In this section we use the following ansatz for nucleon FFs 

)2 \ 2 



F 1 (q 2 ) = F 2 (q 2 )/fi 



Qo 



(30) 



setting the parameter Ql to 0.71 GeV 2 . This form permits us to carry on analytical calcula- 
tions. Note that this presctiption differs from the ones given above, Eq. (1). Let us rewrite 
the expressions Eq. ( 1231) for ^ as follows: 



AM w ^ 4M 
with (here we use dipole approximation for FFs, Eq. [301 ) : 



(31) 



(Qg) 2 



01 

(Q q ) 



•!• = #L^l[(Q fc ) 2 + g 2 (Q,) + g 2 (Q f . 



(*0*(Q*) 



9/J ' 



(OS 



2\2 



(Q 9 ) = <? 2 " Qo, {Qk) = k 2 -Ql 



(32) 



When integrating on the loop momenta, the $i and $2 terms give origin to UV divergences 



contributions containing La = In jgj in the form: 



v4j(s,m) ~ / dxydy 



1 ~ Ji(t) A A 2 3\ ,x , W A A 2 3 



(33) 
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with Dq } do given in Appendix B, and 



G(s,t) = —Trp'^^pa, ■ -^rR llllalv = -F 2 {q 2 ){2t 2 - 6st) - F 1 (q 2 )(2u 2 + 8s 2 + 10tM 2 



F(s,t) = -^TrPil^luPUr, ■ ^TrR'j fl j a [q,jA = -F 2 { 



8Mt 2 + 



8stu 



IQ^MF^q 2 ). (34) 



The $2 and $3 terms contain IR divergences. However, in both regions where IR divergences 
are present i.e., k — > 0, k — > q, the sum of the contributions $2 + $3 converge. As for UV 
contributions, we note that the quantity: 



A 2 



0/ r 



Af(s, u) - Aj(u, s) - 3(u - s)(F 2 + F 1 )(F 1 - 1)L A , L A = In — . 



(35) 



is finite at the limit A — > 00. It results in the replacement A 2 = M 2 , L\ = 0. The explicit 
expression for a,f r in terms of three- fold integrals is given in Appendix B. We note that UV 
divergences associated with $ 2 are canceled due to the symmetry F(s,t) = F(u,t). 
The relevant contribution to the differential cross section can be written as: 

da p dds ot 



dn 



dQ it 



Df, D f 



a fr 

2(B t + B box ) 



(36) 



VI. RESULTS AND DISCUSSION 



The differential cross section with two photon exchange and the relevant soft photon 
emission is given by: 



da da Be 



+ 



dax da 



B 



dn dn dn dn 

Our result for charge asymmetry (T5]) has the form: 

a 



+ 



daF 
~dfl 



A 



odd 



, 1, 2AE _. E 
2 In -In + £>(—, 



D 



tbs 



D 



(37) 



(38) 



where D ths denotes the analytical part 



D tbs = lnxlnp-21n 2 p + Li 2 (l- — J - Li 2 ( 1 - -) 

\ px J \ xJ 



B t 



B t + B box 
+B box lnpln(4r 



1 2 s 1 2 u 2 
m — - — In — - — 71 

M 2 M 2 



2Li 2 1 + 



+ 2Li 2 1 + 



M 



2\ 1 



II 



[1 _ p( s ^ u )] [ A (s, t)ln^ + B(s, t) - t 2 (d 1 F 1 - d 2 F 2 )] } 



(39) 



We note that charge asymmetry A is finite in zero electron mass limit, but contains the 
soft photon emission parameter AE/M. 
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The function D t b s , which can be calculated analytically and the function Df, Eq. (|39|) . for 
which a numerical integration has been performed, have been calculated for several values 
of angles and energies. Both contributions are larger (in absolute value) at large E/M. 

The ansatz which reflects the possibility to separate the QED and the strong inter- 
actions contribution to FFs and that assumes Fqed = 1 is useful as it allows to perform, 
at least partly, analytical calculations. However, it can not be considered exact. A better 
parametrization of QED and strong interaction contributions to FFs is : 

„ . . „ . . „ . . 1 + x 2 „ , . 2x -q 2 Q 2 ( GeV) 

(40) 

Therefore, in order to find the global correction due to two-photon exchange, one should 
weight the individual contributions Dtbs and Df, by the ratio of strong and EM FFs. This 
can be done by multiplying D t b s by the factor Fiq(x). The total contribution is therefore: 

D^V = F 1Q (x)(D tbs + D f ). 

Tabulated values of the numerical results are presnted in Table [H One can see that the 
two photon contribution to the asymmetry is of the order of percent, keeping in mind the 
multiplicative factor oc/tt. The behavior is smooth in the considered kinematical ranges. 
Singularities are expected for 6 — e = 1, due to symmetry properties of the 27 exchange 
fl. 

Taking into account the factor a/n, the corrections do not exceed 1%, in the considered 
kinematical range. 



VII. CONCLUSIONS 



Charge asymmetry in electron proton elastic scattering contains essential information on 
the contribution of 27 exchange to the reaction amplitude. This amplitude can shed light on 
Compton scattering of virtual photons on proton. It contains a part corresponding to proton 
intermediate state, which carries the information on proton FFs. Another term corresponds 
to excited nucleon states and inelastic states such as Ntt, N2tt, NNN. Their theoretical 
investigation is strongly model dependent. 

Our main assumption about the compensation of pure strong interaction induced contri- 
butions to FFs and inelastic channels allows us to avoid additional uncertainty connected 
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with inelastic channels. 

Our choice of photon form factors, (see f. fl30l) ) in nn physical one. The physical case 
corresponds to 

Ge = \ Gm = (QFW (41) 

The aim of the choice is the simplification of the analytical calculation. 

The assumption about the possibility to omit Ai nt in (9) was proved to hold in frame 
of QED, and for ep elastic scattering for the kinematics of almost forward scattering. The 
application to large angle scattering requires, in principle, a rigorous proof. 

The parameterization of the NN*^* vertex is also approximated, since one of the nucleon 
is off mass shell. Nevertheless, they can be estimated to be small for Q 2 < 5 GeV 2 and in- 
cluded in the sources of theoretical uncertainties. The reliability of our assumption can be 
estimated from the ratio of cross sections of pion and nucleon antinucleon pair photopro- 
duction. The uncertainty of our results does not exceed 10%. 

Another interesting question is if such compensation is present for the annihilation chan- 
nel. The measurement of charge asymmetry is, in this case, associated with polar angle odd 
contribution to the differential cross section. 

Similar effects of charge and angular asymmetries can also be due to Z-boson exchange, 
but such contribution is small for moderate-high energy colliders. The ratio of corresponding 
contributions can be evaluated as: ~ {^9v9as) / (aMf) < 5 ■ 1CT 3 for s < 10 GeV 2 (gy (^) 
is the vector (axial) coupling constant of the Z boson to fermion). 

The analytical calculation of 2 / y amplitude with FFs encounters mathematical difficulties. 
In Ref. |8j the results for the box amplitude with arbitrary FFs was investigated. Similar 
attempt was done by A. Ilichev 9j. These works use different approaches to include FFs, 
however the numerical results are close to ours, and show that two photon contribution can 
not be responsible for the discrepancy in recent FFs measurements. 

Other works [lOj devote much attention to the excited intermediate states as A and N* 
resonances, introducing additional uncertainties. In our approach, excited states should not 
be included, as they correspond to poles in the second physical sheet. 

Our numerical results show that charge-odd correlations are of the order of percent, in 
the kinematical region considered here. Such value is expected to be larger at larger q 2 
values and could be measured in very precise experiments, at present facilities. 
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IX. APPENDIX A. LIST OF NECESSARY INTEGRALS. 



We give here a list of scalar, vector and tensor type loop momentum integrals with three 
denominators (k), (e), (p): 

f d kl;k^', kfjky 



I 



m 2 (k)(e)(p) 



Zi] Z 2 pi fl + Z 3 p^; 



Za9w + Z 5 p^p lu + Z 6 pnp u + Z 7 (pi^p u + pivPp). 



(42) 



Standard procedure of joining the denominators leads to integrals of the form 



drr / 2ydy 



1, k^, k^k v 



[(k - yp x ) 2 - y 2 p 2 x - A 2 (l -y)] 3 ' 



(43) 



o o 

with p x = xpi — (1 — x)p,pl = m 2 x 2 + M 2 (l — x) 2 + S\x{l — x), S\ = —s — iO. Further 
integration on loop momentum (see (37)) with A-is the UV cut-off parameter, leads to: 



Zi 

Li s 
Z 2 
Z 4 

z 5 
z 6 
z 7 



1 

2^ 



1 M 2 
L 2 s --L 2 -2Li s + \n—(2L s + L) 



, L = \n — -, L s = ln — — -nr, 



M 2 



Li 2 (1 + )-mln(l + 

s s 

M 2 

L+(l + r^rlL, 



s + M 2 ' 
3 s 



4 La+ 8 4(s + M 2 ) 



1 

2~s 
1 1 



M 2 

-1 ^ + L s (l- 



M 4 



2s + M 2 
1 1 



2s + M 2 



s + M 2 

S -^L S -1 

s + M 2 

M 2 

1 + —L K 



(s + M 2 } 



's + M 2 ^' 



;) + L 



s + M 2 ' 
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Integrals with denominators (k)(e)(p) can be obtained from ones given above by replace- 
ments pi — > —p'i, p — > —p' with the same coefficients. Integrals with denominators (k)(e)(p) 
can be obtained from ones given above by replacements p\ — > —p[, p — > p, and coefficients, 
which can be obtained from the ones mentioned above by replacement s — > w. 
Box type integrals defined as 

k 1; /c^jj k^ky 



Y\, I^iu / ~r 
J %1 



-■ I^Y 2 A^ + Y 3 P^, 



in" (0)(q)(e)(pY 

V = ^ + Y 5 A^A„ + F 6 P M P, + Y 7 (P^ + P,A M ) + F 8 g M Q„ (44) 



with 



Q = 



Pi - Pi 



p = 



Pi + p[ 



A = 



P + p' 



(45) 



Explicit expressions for are: 



V 2 
Y 4 



Y 5 
Y 6 



Y 7 



2, -s - iO , -t 

— In — 7T In t^7> 

st Mm A 2 



1 

2d 



r 



--(F + F Q )-P 2 (F + F A ) 



2d 



T 



--(F + F A )-A 2 (F + F Q ) 



-^(F-G + H P + H A + H Q ) + H A (~ - P 2 ) - H Q {- T - - A 2 ) 



+2Q 2 Y 3 (P 2 + r) + P 2 G A - A 2 G Q + 2Q 2 A 2 Y 2 



1 

Vd 
1 

rd 



+ 2 (p 2 ) 2 (H A - 2g 2 r 3 - g a ) + ip 2 a 2 )(h q + 2Q 2 y 2 - g q ) 

- T -A 2 {G -F-H P -3H A + 6Q 2 Y 3 ) + (P 2 A 2 + ^)(H A - 2Q 2 Y 3 - G A ) 



-A 2 (H Q + 2Q 2 Y 2 -G Q ) 



1 

rd 
1 



T -H - A 2T -{H Q + 2Q 2 F 2 - G Q ) + 2P 2 ^(# A - 2Q 2 r 3 - G A ) 
-t(H a - 2Q 2 Y 3 + H P + \f - X -G) + A 2 (H Q - 2Q 2 Y 3 - Gq) 



Q 2 r 

-P 2 {H A - 2Q 2 Y 3 - G A 



t = 2PA, d = P 2 A 2 -—, H = F-G + H P + 3H A + 6P' Z Y 3 . 



(46) 
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The quantities entering here are: 

1 



M Z F Q 

F A 
F 
Ho 



Vt(1 + t) 



1 , 2 ~t 7T 2 

2 ln ^ + ¥ 



vr 2 + ln(4r) lux + Li 2 (-2y/rx) - Li 2 

1 



2VT 



1 

2^ 



2 In 



-s-iO 



-t 



In 



Mm M 2 
-s-iO 



In — -ln z - 



4M 2 (1 +r) 
■s-iO, 



-tF fl - 2 In 



M 



x 
-t 
M 2 



M 2 

:) + 2 ln^ — + 2Li 2 I 1 + 

M 2 m V s 



M 2 



(47) 



s + M 2 

The explicit expression for QED box-type Born amplitude (see Eq. (21)) is a b (s,u) = 
Ab(s, u) — Ab(u, s) with 

A b (s, u) = (F 1 P q - F 2 Q q )F Q + (FiP A - F 2 Q A )F A + (F 1 P G - F 2 Q G )G Q + 
(F 1 P H - F 2 Qh)Hq + (FiPp - F 2 Qp)F + (FxiV - F 2 Q Y )Y 1 , 



with 



[M 2 ( S - M ) + s(3 S - M )]; g, = -; P A = — - + — ; 



t 2 

MH su 



M 2 t tu _ 

— ; ?r; Qg 



M 2 t st 



tM 2 



St - P 



M 2 ts 



s(s 2 + u 2 



4 8' 

sM 2 sis - u) 

h — -: 

2 4 

st 2 



(48) 



X. APPENDIX B. 

The contributions from the uncrossed box-type Feynman amplitude, Af, can be written 
as the sum of terms associated to $i, $ 2 , $3: 



A/ = A/1 ~~ AM ( F *M A f* + Wo) 2 ^/ 3 ) , a /r = (1 - P(s <- u))A f . 
We follow the procedure described in Appendix C. 

The <3> 2 and $1 terms are calculated using the relation: 



A: 2 -g 2 

(k - q y 



= l-2q t 



(49) 
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in order to eliminate explicitly the divergence at k — > q. For An we have: 



Afi = t / dxydy < G(s,t) 



M 2 3\ | By' 
d 2) d 



with 



and 



tj dxydyz 2 dz \s + 2N{zP ) (j^ + t(2z + yz) 



S = tG(s, t)(2z + yz) + 2(zy - z)[2t{s 2 + tM 2 )Fi - 2st 2 F 2 



N(b) = ^Trp'tf^pt - b)j v pij v ■ ^TrR-f^(p + b + M)j v 



and A, B given above (Eq. (1321)). The expression for Af 2 is: 

Af2 = 

+ 



dxdyyz 2 dz ^ F(s, t) 
1 

Do 



ln^-- 3 -)- q 2 yL 1 + q 2 (Q q ).y 2 L 2 



q 2 yL 2 - 2q 2 (Q q )fL 3 



[M(yp x ) + M(yp' x )} 



2 / dxdyyz 2 dz 



1 q 2 z 



7T + V^2 - q\Q q )rJi 



t[2z + zy] 



2a 2 

~„2 <T l C„2 



q a [F a (zP )+F a (zP^ 



--2zq 2 J 3 + 6q 2 (Q q )z 2 J 4 



[M(zP ) + M(zP^)] 



with Ji, Li given in Appendix C and 



M ( b ) = ^Trp'il^Pi - b)j v pijr, ■ -TrRypip + b + M)[q, 
M (6) = -Trp'^^ + b) lu p llv ■ ~TrR[q, 7 „] (p'-b + M ) 7 „. 

1 1 

F a (b) = ^Trp 1 7 M 7 a 7 l/ pi7„ • -TrRry^ + b + M)[q, 7„] - 

|Trpi 7m (pi - 6)7^i7r? • jrr^algj,] + (g - b) a F(s,t); 



and similar expression for F a 
The form of A j 3 is: 



A /3 = / rfx^(^) 2 ^{-[^( 2 Po) + ^(^o)]J3 + 6[G( 2 P ) + G(^P ( ;)]J 4 } 



18 



with 



G(k) = ^Trpi 7 M (p! - k)^ v pi^ v ■ ^TrR-f^p + k + M)[q-k, j v ]; 

G(k) = ^Trft-y^ + kh^ ■ -TrR[q - k, 7/ J (p' - fc + Af ) Tl/ ; (54) 

#(&) = ^ Tr P'i7^7 CT 7^i7r, | ~TrRj^ a [q - 6, 7 „] + ^TrR-f^p + b + M)[j a , j v ] J 
1 - 1 

#(&) = ^ Tr P'i7^7<77^i7r, j-^rRfy^tp' - & + M )7^ - - 6, 7i/]7<r7* j 



+|Trp' 1 7 M (pi / + b)-f u pa v ■ - Tr #[7^] 7^7*,. 



(55) 



Infrared singularities (divergences of Af2, Af 3 at y — > 0, z — > 1 are mutually compensated 
in the sum Af. 

XI. APPENDIX C. 

Let us describe the procedure employed for compacting the denominators and for the 
loop momentum integration. Taking T>(z) = az + bz, z = 1 — z, the Feynman prescription 
for compacting the denominators leads to: 

1 f l dz 1 f 1 2zdz 1 f 1 6zzdz 



(56) 



ab J Viz) 21 a?b J Viz) 3 ' a?b 2 J V{zf 

1 f 1 3z 2 dz 1 f 1 \2z 2 zdz 1 f 1 30(zz) 2 dz 



a 3 b J V{z) v a 3 b 2 Jo V(z) 5 ' a 3 b 3 J V(z) 6 
Applying (|56|) to the set of denominators: 

(e) = k 2 - 2Jfepi, (e) = k 2 + 2kp[, (p) = A; 2 + 2£;p; 

(p) = k 2 - V*, (*) = K 2 , (k) = (q~ k) 2 , (Q k ) = k 2 - Ql 

one obtains: 

1 r i 2dxydy 1 f 1 2dxydy 



f 1 
Jo 



(epfc) J0 [fc-y Pa: )2-P ]3' ( epk ) Jo [(k-P ) 2 -Vo} 3) 
1 „! 2dxydy 1 f 1 2dxydy 

JO 17 1. DA 2 T) 13' l^/O^ ~~ / f/' JU ^2 .7 13' v ' 



(epfc) J0 [(A; - P^) 2 - V0} 3 ' (epQ fc ) 7 [(fc - yp,) 2 - rf 
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with £> = y 2 p 2 x ; d = y 2 p 2 x +yQo, Po = yPx+yq, Po = yp' x +yq, p x = xpi-xp, p' x = xp '-xp\. 

The following relation holds: P 2 = P^ = y 2 p\ + yq 2 , pi = p' 2 = x 2 M 2 — sxx. 

We do not develop further (explicitly) the denominators which contain (e)(p), because 
the corresponding results can be obtained by those depending on (e)(p) under replacement 

Px -> Px'- 

In a similar way one obtains: 

1 f 1 Qyydxdy 1 f 1 Qyydxdydt +t-r> 2 

(epQ 2 ) J [(k - yp x ) 2 - d o r (ke P Q k ) J [(k - yp x ) 2 - dtf' 1 ° + yWo ' 

1 f 1 24dxtdtyy 2 dy 1 f l 6ydydxz 2 dz _ 2 

WpQl) = J [(k - yp x ) 2 - d t f' WQk) = Jo [(k-zP ) 2 -V i r i = a + ^o; 

1 f 1 Qydydxz 2 dz 2 2 _ 2 

t — ttt = / 777 75 — u ;a = z p o - z yq ; 

1 Z" 1 24ydyz 2 zdzdtdx _ 2 

^pkkQ^) = Jo [{k-zPo) 2 -V 2 r V2 = a + Z ^ 

1 Z" 1 120ydydxtdt(zz) 2 dz 

(epkkQ 2 ) = Jo [(k-zPo) 2 -V 2 f ■ ( } 

After replacing the loop momenta by fc — 6 — > /c, a symmetrical integration on k is performed. 
For polynomials in N(k) of order in k < 3 one finds: N(k + b) = \K 2 (b) + AT(6). Applying 
Wick rotation leads to: 

r K 2 dn . A 2 3 f dn 1 
i 77^-773 =ln ^- o 



n 2 -df d 2' y (K 2 -d) 3 2d' 

r K 2 (ift; 1 f dn 1 

J (k 2 -dy = ~m ] J (k> - dy = od 2 ' 

K 2 dK 1 f dn 1 



„ K diK 1 f 

J (/t 2 - rf) 5 = i^ 2 ' y 



(/t 2 -d) 5 12rf 3 ' 



f « 2 (ifi; 1 f dn 1 

J {k 2 - dy = ~sod? ; y (k 2 - d) 6 = 2o7i4 ; ( ' )9j 



with 



d K 2 i ■) ' } 2 ° 2 2 

UK = ~ ^ = /tgfl/tg, /t = Kg — /tg -|- -\~ AC 2 -|- Kg, 

and K e - euclidean four vector. Note that the integration of the Feynman parameter t can 
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be provided explicitly, as it enters only in phase volume: 

11 dt 1 , d T f 1 dt 
m — ; L 2 



tdt 



1 dt 

o U 'i 
1 tdt 

O U 2 



(VQo) 2 

a + Vx 
2a 2 V\ ' 

j J\ 



In 



td 



d 2 t V d ' 
d 



2aVf 



2 ' 



1 tdt _V 1 + 2a 



(60) 
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E/M-9 


30 


60 


90 


120 


150 


1 


1.16 


0.58 


-0.72 


-1.29 


-1.81 


2 


3.71 


1.19 


-0.22 


-0.96 


-1.61 


3 


4.23 


1.54 


-0.15 


-0.97 


-1.73 


4 


4.46 


1.86 


0.32 


-0.62 


-1.35 


5 


4.52 


2.43 


1.12 


-0.13 


-0.99 



TABLE I: Numerical values of I? as a function of E/M and 6, with dipole parametrization of the 
form factors. 
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